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The Peierfs instabifity in multi-channei metai nanowires is investigated. Hyperscaiing reiations 
are established for the finite-size-, temperature-, and wavevector-scaiing of the eiectronic free energy. 
It is shown that the softening of surface modes at wavevector q = 2fc F , v leads to critical fluctuations 
of the wire's radius at zero temperature, where fc F , v is the Fermi wavevector of the highest occupied 
channel. This Peierls charge density wave emerges as the system size becomes comparable to the 
channel correlation length. Although the Peierls instability is weak in metal nanowires, in the sense 
that the correlation length is exponentially long, we predict that nanowires fabricated by current 
techniques can be driven into the charge-density-wave regime under strain. 
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I. INTRODUCTION 

Already long ago, Frohlich 1 and Peierls 2 pointed out 
that a one-dimensional metal coupled to the underly- 
ing lattice is not stable at low temperatures. Electron- 
phonon interactions lead to a novel type of ground state 
with a charge-density wave (CDW) of wavevector 2k p 
(see Ref. 3 for a review). This state is characterized 
by a gap in the single-particle excitation spectrum, and 
by a collective mode with an associated charge density 
~ po + pi cos(2fcp2:), where po is the unperturbed elec- 
tron density of the metal. Of particular interest are in- 
commensurate systems, where the period of the CDW 
is not simply related to that of the unperturbed atomic 
structure. In that case, no long-range order is expected 
even at zero temperature due to quantum fluctuations. 

In contrast to the usual Peierls systems, metal- 
lic nanowires 4 are open systems with several inequiv- 
alent channels, for which the theory has not yet 
been developed. Interest in the Peierls transition 
in metal nanowires has been stimulated by recent 
experiments 5,6,7 ' 8 ' 9 on nanowire arrays on stepped sur- 
faces. In these systems, interactions between nanowires, 
mediated by the substrate, render the system quasi- 
two-dimensional at low temperatures, similar to the 
dimensional crossovers commonly observed in highly- 
anisotropic organic conductors. 3 However, individual 
free-standing metal nanowires 4 represent true (quasi)onc- 
dimensional systems, in which the intrinsic behavior of 
the Peierls-CDW can be studied. 

Due to quantization perpendicular to the wire axis, 
electron states in metal nanowires are divided into dis- 
tinct channels, which are only weakly coupled. Each 
channel has a quadratic dispersion relation, and starts to 
contribute at a certain threshold energy, i.e., the eigenen- 
ergy E n of the corresponding transverse mode. This re- 
sults in a sequence of quasi-one-dimensional systems with 
different Fermi wavevectors 



k F , n = ^-^{E F -E n ). (1) 



The channel Fermi wavevectors are generically not com- 
mensurate with the underlying atomic structure in 
nanowires with more than one open channel. 

The Peierls instability is weak in metal nanowires, 10 
so that the system is close to the quantum critical point 
at which the transition from Fermi liquid behavior to a 
CDW state occurs. In the vicinity of the quantum critical 
point, the system exhibits an additional length scale, the 
correlation length 

Here, the greek index v labels the highest occupied chan- 
nel, in which an energy gap 2A„ opens and v F , „ = 
(h/m e )k Fjl/ is the Fermi velocity of this subband. Within 
a hyperscaiing ansatz, 11 the singular part of the energy 
is expected to scale like E s [ ng /L ~ £,J 1_Z , where L is the 
wire length, and the dynamic critical exponent takes the 
value z — 1. Thus, E s i ng /L ~ Al/hv F . v . Near the sin- 
gular point, we indeed find that the electronic energy is 
given by 

3** « *L Y (z vSq ,k&), (3 ) 
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where Y(x,y,z) m ln(max{x, y, z}) is a universal and 
dimensionless scaling function, 5q = (q — 2fc P) „) is the 
detuning of the perturbation wavevector from its critical 
value 2k^ Vl Lt — %v? iV /kBT is the thermal length at 
temperature T, and n v = 1 or 2 is the degeneracy of the 
highest open channel v. This universal scaling behavior 
is quite different from that of a closed system with pe- 
riodic boundary conditions, 12,13 where radically different 
behavior was found for odd or even numbers of fermions. 
Nonetheless, the correlation length £ was also found to 
control the finite-size scaling of the Peierls transition in 
mesoscopic rings. 13 

In this article, the quantum and thermal fluctuations 
of the nanowire surface are calculated in a continuum 
model, 14 where the ionic background is treated as an in- 
compressible, irrotational fluid. In contrast to the semi- 
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classical theory of Ref. 14, which exhibited critical sur- 
face fluctuations only at finite temperature — due to the 
classical Rayleigh instability — the present fully quantum- 
mechanical theory exhibits critical zero-temperature sur- 
face fluctuations at wavevector q = 2fc Fi „. In a finite 
system, these CDW correlations are found to grow in am- 
plitude as the wire length L approaches a critical length 
L c of order the correlation length £„. A similar scal- 
ing is observed as the temperature is lowered, so that 
Lt exceeds £„. Although £„ is typically very large for 
fully-equilibrated structures, consistent with the fact that 
Peierls distortions have not yet been observed in multi- 
channel nanowires, it is predicted that nanowires of di- 
mensions currently produced in the laboratory can be 
driven into the CDW regime by applying strain. 

This paper is organized as follows. Sec. II summarizes 
the (standard) Peierls theory for a one-dimensional metal 
with a half-filled band, and extends it to multi-channel 
wires. A description of the deformation of a nanowire 
through surface phonons is given in Sec. III. The corre- 
lation length £ is introduced in Sec. IV, and the scaling 
relation (3) is established. Sec. V examines the critical 
surface fluctuations, and consequences for different ma- 
terials are discussed. Finally, a summary and discussion 
are given in Sec. VI. 



II. THE PEIERLS INSTABILITY 

Consider the ground state of a one-dimensional lin- 
ear chain of atoms, with lattice constant a and periodic 
boundary conditions: In the presence of electron-phonon 
interactions, it is energetically favorable to introduce a 
periodic lattice distortion with period A = ir/k F ; this 
effect is known as the Peierls instability. 2 The lattice dis- 
tortion opens up an energy gap 2A in the electronic dis- 
persion relation at the Fermi level E F , so that the total 
electronic energy is lowered. If A <C E F , then the gain 
in energy is given by 2 

L Tr2E F /k F [ & \4E F J 2\ \E F ) K ' 

The size of the gap can be extracted from an energy 
balance: Let the increase of the elastic energy due to 
the deformation be given by E cost /L = ab 2 , where b is 
the amplitude of the distortion. On the other hand, A is 
defined through the matrix element of the perturbation 
coupling states with longitudinal wave vector ±k F ; it is 
linear in b and we can set A = A b. By finding the 
minimum of 5E(b) = E ga i n + E cost , we can derive the 
optimal value of the distortion amplitude b and from this 
derive the size A of the gap, 

A = AE F vip(-2iraE F /A 2 k F ) . (5) 

The analysis of the Peierls instability in a one- 
dimensional metal can be extended to the case 



of a multi-channel system. While other instabili- 
ties of the Fermi liquid — induced by electron-electron 
interactions — compete with the Peierls instability in 
purely one-dimensional systems, 3 their importance de- 
creases as the number of channels increases, 15 so that 
electron-electron interactions can in a first approxima- 
tion be neglected in multi-channel nanowires. Moreover, 
electron-electron interactions are strongly screened in s- 
orbital metal nanowires with three or more conducting 
channels. 16,17 Hence, including electron-electron inter- 
actions in the calculation will not lead to any qualita- 
tive change (except in the limit of very few conduction 
channels) whereas electron-phonon interactions do. 18 We 
therefore consider only electron-phonon coupling in this 
article. 

For any given channel n, a perturbation of wavevec- 
tor q = 2k^ n will open a gap 2A„ at the Fermi sur- 
face in the energy dispersion of this channel. The energy 
gain -Egain.n and gap size 2A„ are then given by Eqs. 
(4) and (5), respectively, with k F replaced by fc Fj „ and 
E F replaced by h 2 k 2 n /2m e . The greatest effect will be 
seen close to the opening of the highest occupied channel 
(i.e. n = v), where the channel Fermi wavevector fc F „ is 
small. Note that the same perturbation (with q = 2£; Fj! ,) 
will also modify the dispersion relations of lower-energy 
channels n' with E n i < E v , but due to the finite spac- 
ing of the threshold energies, this modification will occur 
within the Fermi sea and there will be little net effect. 10 

The standard Peierls theory uses periodic boundary 
conditions for the wave functions. Its extension to the 
multi-channel case cannot be directly applied to a metal- 
lic nanowire of finite length L. The nanowire is part of a 
much larger system including the leads, and the longitu- 
dinal wavevectors k n in the subbands are not restricted 
to multiples of as in the case of an isolated system 
with periodic boundary conditions. Therefore a pertur- 
bation of the nanowire with wavevector q does not only 
couple states k n and k' n that exactly obey k n = k' n + q. 
Instead the state k n is coupled to a range of fc^-states 
proportional to l/L. The dispersion relation remains 
smooth while, with increasing wire length, it develops 
a smeared-out quasi-gap, and only in the limit L — > oo 
do we recover the Peierls result with a jump at fc Fj „. 10 

III. SURFACE PHONONS 

We describe the wire in terms of the Nanoscale Free 
Electron Model (NFEM), 19,20 treating the electrons as a 
Fermi gas confined within the wire by a hard-wall poten- 
tial. The ionic structure is replaced by a uniform (jel- 
lium) background of positive charge, and we assume that 
this ionic medium is irrotational and incompressible. The 
approximations 14,16,20 of the NFEM require strong der- 
ealization of the valence electrons, good charge screening, 
and a spherical Fermi surface, conditions met in alkali 
metals and, to a lesser extent, noble metals such as gold. 

In this continuum model, the ionic degrees of freedom 



are completely determined by the surface coordinates of 
the wire. Let us consider an initially uniform wire of ra- 
dius i?o and length L which is axisymmetrically distorted. 
Its surface is given by the radius function 



R(z,t) = Roll 



(6) 



where the time-dependent perturbation is written as a 
Fourier series with coefficients b q (t). Since R(z, t) is real, 
we have b q = b*_ q and we require that the volume of the 
wire is unchanged by the deformation. Other physically 
reasonable constraints arc possible and will be discussed 
in detail in Sec. V. 

The kinetic energy of the ionic medium is given by 14 



E, 



kin 



= L^2m(q,R a ) 

q>0 



R 



db q (t) 
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(7) 



Here, the mode inertia m is a function of wire radius and 
phonon wavevector, and reads 

27ri? / (<7i?o) 



m(q,Ro) = pi, 



qh(qRo) 



(8) 



where pi on is the ionic mass density, and Iq and 1\ are 
the modified Bessel functions of zeroth and first order, 
respectively. 21 Considered as a function of q, the mode 
inertia has a singularity ~ 1/q 2 at q = 0, and is mono- 
tonically decreasing, with m ~ 1/q for large q. 

Within the Born-Oppenheimer approximation, the po- 
tential energy of the ions is given by the grand canonical 
potential SI of the confined electron gas. A linear stabil- 
ity analysis of cylindrical wires determining the leading- 
order change in SI due to a small perturbation was re- 
cently presented in Ref. 10: (SSI is quadratic in the Fourier 
coefficients b q of the deformation, and can be written as 



5CI 

T 



J2\b q \ 2 a(q,R a ,L,T) 

q>0 



(9) 



where the terms of order Xf / L include nondiagonal con- 
tributions which can be neglected if the wire is long 
enough. 22 The explicit analytical expression for the mode 
stiffness a is given in App. A. Here, we are interested in 
the general behavior of the mode stiffness as a function 
of the perturbation wavevector q for a given radius Rq, 
length L, and temperature T: a{q) is a smoothly increas- 
ing function, with a L- and T-dependent dip at q = 2fc F , „, 
where v is the index of the highest open channel (see App. 
A). Therefore, we formally split a into a smoothly vary- 
ing term and one containing the singular contribution of 
channel 



& — "smooth 



» 



(10) 



The smooth part can be thought of as the sum of a Weyl 
contribution, 14 which describes the effects of surface ten- 
sion o~ s and curvature energy 7 S , plus an electron-shell 
correction 14 

"smooth = - 27r.Ro °~s + 2-kRI (er s R q - j s ) q 2 + a shcn , (11) 
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FIG. 1: (a) Surface phonon dispersion relation oj(q) for a wire 
of length L — L c , where L c is the critical length (see text). 
Different wire lengths are compared in the inset showing the 
vicinity of q — 2fc Fj „. The curves are offset horizontally by 
0.05 units for clarity, (b) Temperature dependence of the 
minimum in u>(q) for a wire of length L = L c (curves are offset 
horizontally by 0.1 units). The wire radius is J?o = 4.5fc^ for 



all curves, thus L c 



6500fc" 



whereas the singular part describes the onset of the 
Peierls instability in channel v. At zero temperature, 
it is given by 10 



2k FU +q 



In 



2fc F 



-9 
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h 2 nq 

-F((2k F>u +q)L 



(12) 



F(|2fc Fj v —q\L 



where F(x) = Ci(x) — sin(a;)/x and Qi{x) is the cosine 
integral function. 23 The finite temperature mode stiffness 
is evaluated numerically by computing the integral 



a(T) 



dE 



(13) 



where f(E) is the Fermi function and a(E) is obtained 
from the zero temperature results by replacing fc F: v by 
k E , v = [^p^(£ — E v )]~ 1 / 2 . 

Combining the kinetic energy (7) and potential en- 
ergy (9) yields a Hamiltonian for surface phonons, H ph = 
2 ), with frequencies 



J2 q nuj q {a\a q + 1 



w(q,R ,L,T) 



l a{q,Ro,L,T) 
R%m(q,R ) 



(14) 



These axisymmetric surface phonons correspond to the 
longitudinal acoustic mode of the nanowire. From Eqs. 
(8), (10) and (14), we infer that o->(q) is a smoothly in- 
creasing function of q, with a dip at q = 2k FjU . As 
expected for acoustic phonons, it is linear for small q. 
The i-dependent softening of the phonon modes with 
wavevector 2fc F v defines a critical length L c for which 
uj(2k F , v ) = 0. 

A plot of Lu(q) for L = L c at zero temperature is shown 
in Fig. 1 (a). The inset shows a close-up of the minimum 
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and compares different wire lengths, where the curves 
are horizontally offset for clarity. The temperature de- 
pendence of u> is illustrated in Fig. 1(b), concentrating 
on the vicinity of q — 2fc Fj v (again, the curves are hori- 
zontally offset). With increasing temperature, the dip in 
lu disappears, and the curve becomes smoother. 



IV. SCALING RELATIONS 

The opening of the Peierls gap in subband v introduces 
a new energy scale, given by the gap 2A„ for an infinitely 
long wire. A„ is determined by the matrix element of 
the perturbation coupling states with longitudinal wave 
vector ±fc Fi „ and is linear in the distortion amplitude b q . 
The perturbation potential matrix element is calculated 
in App. B, and we find that A„ = 2E u b q . 

On the other hand, the energy cost for creating a sur- 
face modulation with wavevector q is determined by the 
smooth part of the mode stiffness, so that SE cost /L = 
&smooth(q)\b q \ 2 ■ Following the arguments of Sec. II we 
can now calculate the length scale £„, obtaining 



liA7 



1 



4fc F 



exp 



2m. 



Trkv i/C^smooth 



2k v EI 



(15) 



where we have used Eq. (5). Introducing the cor- 
relation length £„ allows us to derive the finite-size-, 
temperature-, and wavevector-scaling of the electronic 
free energy near the critical point of the Peierls insta- 
bility, given by q = 2/c Fj „, L — > oo, T = 0, and weak 
electron-phonon coupling. 

Finite-size scaling — First we examine the mode stiff- 
ness at zero temperature and for q = 2fc Fj „ as a function 
of wire length. Starting from Eqs. (10) and (12), we take 
the limit q — > 2fc Fi „ and get 



a(L) 



T = 



a SBl ooth - ^ K " Ev [ln(4fc P „L) - ci] , (16) 
irnv F v 



where Ci = l-j E + F(4fc F , V L) » 0.42. Here ~/ E ^ 0.577 
is the Euler-Mascheroni constant, and we have used the 
fact that F(x) <C 1 for x ^> 1. This expression for a(L) 
can further be simplified by the use of Eq. (15), so that 



a(L) 



q=2k F) t 
T = 



7rhv v u 



„4 



ci 



(17) 



This defines the critical length L c for which a takes the 
value zero, 



6 « 1-52^ 



(18) 



Note that the critical length is of the same order of mag- 
nitude as the correlation length. 

Wavevector scaling — Now let 8q = (q — 2fc F „) be the 
detuning of the perturbation wavevector from its critical 
value 2fc Fi „. At zero temperature and in the limit L — > oo, 




y/z = l t /l 



FIG. 2: (color online) Cross-over from T-scaling to L-scaling: 
The scaling function Y(a:, y, z) is plotted as a function of the 
ratio y/z = Lt/L for fixed values of y — (, U /L (red curves) 
and 2 = (,v/Lt (green curves) at x = Sq^ = 0. For all curves 
Ro = 8fcp\ 



we expand the mode stiffness as a function of Sq and 
obtain 



a(Sq) 



L — L c 
T=0 



HtJq\ + 0(Sq) . (19) 



Again, the result was written in a compact form by the 
use of Eq. (15). 

Temperature scaling — Finally, we examine the effect of 
finite temperature for a wire of infinite length and a per- 
turbation of wavevector 2fc Fi „. The main effect of finite 
temperature is to smear out the Fermi surface, so that the 
critical wavevector q is detuned by Sq E = 2(fc FiI , — fc Ejl ,) 
with k E , u = (E - Ey)]- 1 ' 2 . Starting from the scal- 
ing behavior for finite Jg-detuning, Eq. (19), we cal- 
culate a(T) for small T from Eq. (13) by linearizing 
E — Ep w Tiv FjU (k Eyl/ — fc F „) and find 



a(T) 



L = L C 
q=2kp 



Trhv F „ 



In 



+ c 2 



. (20) 



with a numerical constant ci = — 7_e + ln7r w 0.57. 

Combining the three scaling relations for length, per- 
turbation wavevector, and temperature, Eqs. (17), (19), 
and (20), respectively, we prove Eq. (3) for the singular 
part of the electronic energy. 24 

Finite temperature introduces a thermal length Lt — 
/3hv Fj u . Changing Lt has an equivalent influence on the 
system properties as changing the length of the wire L. 
Near the critical point, we observe a crossover between 
finite-T scaling and finite-L scaling, depending on the 
ratio Lt/L (see Fig. 2). The smaller of the two lengths 
dominates the behavior of the singular part of the elec- 
tronic energy. As long as Lt < L, a change of L, even by 
orders of magnitude, has only a small effect on Y, whereas 
the system is sensitive to small changes in Lt and shows 
T-scaling. The situation is reversed for Lt > L: In 
this change in temperature results in only small 

changes of Y, whereas the singular part of the energy 
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FIG. 3: (color online) CDW correlations for various values 
of L/L c at T = for a wire with R = 4.42fc j ; 1 . The critical 
length is L c ~ 1560/c^ 1 . Curves offset vertically for clarity. 



depends strongly on L and shows finite-length scaling. 
These two different cases are illustrated in Fig. 2 by the 
two sets of curves for different fixed values of L and Lt, 
respectively. 

So far we have considered ideal nanowires without dis- 
order. Disordered structures exhibit an additional length 
scale, the electron elastic mean free path I. The scaling 
theory we have derived allows us to predict that the ef- 
fect of disorder is to cut off the logarithmic scaling of the 
Peierls-CDW instability, exactly like the thermal length 
or wire length. We thus infer that the length-dependent 
criterion for the emergence of the Peierls-CDW in metal 
nanowires is given by 



L r < L. L n 



(21) 



where the critical length L c (Eq. 18) is of order the cor- 
relation length . Increasing temperature leads to a de- 
creasing thermal length, and therefore destroys the phe- 
nomenon at sufficiently high T. Increasing disorder has 
a corresponding effect. 



V. SURFACE FLUCTUATIONS 

The softening of the surface phonon modes with 
wavevector q — 2fc P) „ leads to critical surface fluctua- 
tions. Given the mode stiffness a and phonon frequency 
u)(q), the fluctuations about the cylindrical shape are 
given by 25 



TOfl(o)) 

Rl 



-s' * 



hu>(q) 
a(q) 



1 



(22) 

Figure 3 shows the correlations for different wire lengths 
at zero temperature. Fluctuations with q = 2fc F: „ in- 
crease with increasing wire length. Note that these CDW 
correlations may be pinned by disorder, or at the wire 
ends. 3 The correlations shown in Fig. 3 are representa- 
tive of regions far from an impurity or wire end. 



For large z, we can use a saddle point approximation 
to estimate the integral in Eq. (22), and find 



(R(z)R(O)) 



oc cos(2fej„js) Jf (y i21 °g(W£) 1 



(23) 

where Kq is the modified Bessel function of the second 
kind 26 of order 0. Its asymptotic behavior is given by 27 



K (x) 



-7_e - logx/2 
~-e~ x 



forx«l, 
forx>l. 1 ' 



Depending on the ratio of the wire length L to the critical 
length L c , we can distinguish two regimes: for L -C L c , 
the prefactor of z /L in the argument of Kq in Eq. (23) 
is large and the correlations decay exponentially for suf- 
ficiently large z. On the other hand, if L ~ L c , our 
theory predicts a logarithmic decay of the correlations. 
Note that the ratio L/L c determines the crossover from a 
regime where the harmonic approximation about a Fermi 
liquid is valid (L < L c ) to that of a fully-developed CDW 
(L > L c ). The harmonic approximation, which we have 
used in our calculation, breaks down at L = L c , where 
the wire can no longer be treated as a cylinder with small 
perturbations. 

The correlation length £ is a material-specific quan- 
tity, since it depends exponentially on the smooth con- 
tribution (11) to the mode stiffness [cf. Eq. (15)]. As 
discussed in detail in Ref. 28, the material-specific sur- 
face tension and curvature energy can be included in the 
NFEM through a generalized constraint on the allowed 
deformations of the wire 



T)s kp S 



rjc k F C 



const. 



(25) 



Here V is the volume of the wire, S its surface area, and 
C its integrated mean curvature. The constraint dfif = 
restricts the number of independent Fourier coefficients 
in Eq. (6), and allows bo to be expressed in terms of the 
other Fourier coefficients. This results in a modification 
of the smooth part of the mode stiffness (see App. A). 
Through an appropriate choice of the dimensionless pa- 
rameters r] s and rj Cl the surface tension and curvature 
energy can thus be set to the appropriate values for any 
given material. 28 

The upper panel of Fig. 4 shows the correlation length 
calculated using the material parameters 29 for Na and 
Au. For clarity, the plot is restricted to the cylindri- 
cal wires of so-called magic radii that prove to be re- 
markably stable even when allowing symmetry-breaking 
deformations. 28,30 These wires have conductance values 
G/G = 1,3,6, 12, 17, 23, 34, 42, 51, where G = 2e 2 /h 
is the conductance quantum. 31 Arrows point to the posi- 
tions of the most stable wires, defined as those with the 
longest estimated lifetimes, 32 which are near the minima 
of the shell potential for straight cylindrical wires, shown 
in the lower panel of Fig. 4. 

The maximum length of free-standing metallic 
nanowires observed so far in experiments is that of gold 
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FIG. 4: (color online) Upper 
panel: Correlation length f„ 
for Na and Au. For clar- 
ity, the plot is restricted 
to wires of so-called magic 
radii, i.e. wires of conduc- 
tance G/Go = 1,3,6,12,..., 
that were shown to be re- 
markably stable even when 
allowing symmetry breaking 
deformations. 10,30 The lower 
panel shows the shell poten- 
tial, and arrows mark the ge- 
ometries with the longest es- 
timated lifetime. 32 



wires produced by electron beam irradiation of thin gold 
films in ultra high vacuum, 33 for which lengths of L ~ 3- 
15 nm are reported. In those experiments, no sign of 
an onset of the Peierls-CDW was seen, presumably in- 
dicating that L < £. However, since the correlation 
length depends exponentially on the wire radius, we pre- 
dict that nanowires of currently available dimensions can 
be driven into the CDW regime by applying strain. A 
tensile force of order 1 nN can change £ by orders of 
magnitude, and thereby drive the system into the CDW 
regime as soon as the condition (21) is met. Note that 
since ftp is of order 1 A -1 for Na and Au, the ther- 
mal length Lt at room temperature is of order 100 A 
for Na and twice as large for gold, comparable to the 
lengths of the longest free-standing wires currently pro- 
duced. By contrast, the elastic mean-free path can be 
several times as long 4 and electron microscope images of 
gold nanowires 33 show perfectly regular and disorder-free 
atomic arrangements. Surface roughness does not play a 
role in the parameter regime we consider, where electron 
shell effects dominate over ionic ordering. Thus CDW be- 
havior should be observable at room temperature in free- 
standing metal nanowires under strain, in contrast to the 
behavior of quasi-one-dimcnsional organic conductors, 3 
where CDW behavior is observed only at cryogenic tem- 
peratures. 



VI. SUMMARY AND DISCUSSION 



is material dependent. We predict that the Peierls-CDW 
should be observable at room temperature in currently 
available metal nanowires under an applied strain. 



The critical length is shortest in materials whose sur- 
face tension is small in natural units (i.e., in units of 
Epkp). Notable in this respect 29 is Al with a s — 
0.0018Epkp, some five times smaller than the value for 
Au. Although Al is a multivalent metal, it has a very free- 
electron-like band structure in an extended-zone scheme, 
and thus may be treated within the NFEM, although 
the continuum approximation is more severe. We thus 
predict that Al should be an ideal candidate for the ob- 
servation of the Peierls-CDW. 



Our findings on the finite-size scaling of the Peierls- 
CDW in metal nanowires are in contrast with previous 
theoretical studies 12,13 of mesoscopic rings: For spinlcss 
fcrmions in a one-dimensional ring, the Peierls transi- 
tion was found to be suppressed for small systems when 
the number of fermions is odd, but enhanced when the 
number is even. Obviously, no such parity effect occurs 
in an open system, such as a metal nanowire suspended 
between two metal electrodes; we find that the CDW is 
always suppressed in nanowires with L <C £ u . Nonethe- 
less, the correlation length £ was also found to control the 
(very different) finite-size scaling in mesoscopic rings. 13 



In conclusion, we have presented a scaling theory of 
the Peierls-CDW in multi-channel metal nanowires. Near 
the critical point, scaling relations for the L-, q- and In- 
dependence of the singular part of the free energy, which 
drives the Peierls instability, were established. A hy- 
perscaling ansatz was verified and the universal scaling 
function was analyzed, which was found to be logarith- 
mic. The crossover from a regime where the harmonic 
approximation about a Fermi liquid is valid [L < L c ) 
to that of a fully-developed CDW (L > L c ) occurs at a 
critical length of order the correlation length £„, which 



The finite-size scaling of the Peierls-CDW in metal 
nanowires is similar to that of the metal-insulator tran- 
sition in the one-dimensional Hubbard model. 34 In both 
cases, there is no phase transition in an infinite system, 
because the critical electron-phonon coupling and on-site 
electron-electron repulsion is + in each case. However, 
for parameters such that the gap 2A in an infinite system 
is sufficiently small, the system is close to a quantum crit- 
ical point, 11,34 and the crossover from Fermi- liquid be- 
havior to a Peierls-CDW or Mott insulator, respectively, 
can be described within hyperscaling theory. 
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APPENDIX A: LINEAR STABILITY ANALYSIS 



where a factor of 2 for spin degeneracy has been included. 
Finally, the grand canonical potential is related to the 
density of states D(E) by 



n 



k B T / dE D(E) In 



(A2) 



This appendix gives details on the linear stability anal- 
ysis (cf. Rcf. 10) which determines the leading-order 
change in the grand canonical potential £1 of the elec- 
tron gas due to a small deformation of a cylindrical 
nanowire. Since the nanowire is an open system con- 
nected to macroscopic metallic electrodes at each end, 
it is naturally described within a scattering matrix ap- 
proach. The Schrodingcr equation can be expanded as a 
series in the perturbation, and we solve for the energy- 
dependent scattering matrix S(E) up to second order. 
The electronic density of states can then be calculated 
from 

D(E) = ^{sHE)^- d ^S(E)}, (Al) 



where ks is the Boltzmann constant, T is the temper- 
ature, and [i is the chemical potential specified by the 
macroscopic electrodes. We find that the change Sfl due 
to the deformation of an initially axisymmetric geometry 
in leading order is quadratic in the Fourier coefficients b q 
of the deformation (cf. Eq. 6) and can be written as stated 
in Eq. (9), defining the mode stiffness a(q, Ro, L, T). It 
is convenient to decompose a into three contributions, 
ol = oyiag a nond + <^con- The coupling between channels 
mediated by the surface phonons determines ayiag, com- 
ing from scattering into the same channel (Eq. A3), and 
anond, coming from scattering between different channels 
(Eq. A4): 



12E n k F n 



a diag {q,R ,L) = -V ]e(E F -E n ) 

IT L — ' 

n 

a n ond(<7, Ro, L) = fn,n>Q(EF — E n ) 

7T ' 4 



In 



q 



2&f, n Q 



16fc F 



E n E n i AE n E n i 

H In 



F{{2k Fn + q)L)+F{\2k F , n -q\L) 
q 2 + E n , -E n + 2qk Fi n 



(A3) 



' TP J? 



+ 



q 2 + E n , -E n - 2qk F , 
+e(E F -E n ,) iEnEn ' (F(\q - fc P; „ - k F . n ,\L) - F((q + k F , n + k Fn ,)L)) 



a con {q,R ) = - > <d{E F -E n )AE n k F , n — . 

^* 1 - r) s /(Rok F ) 



(A4) 
(A5) 



Here /„,„' = 1 for two channels having the same az- 
imuthal symmetry and /«.„' = otherwise. The function 
F(x) = Ci(a;)— sin(x)/x smoothens the logarithmic diver- 
gences so that «dia g and a non d are continuous functions 
having minima of length-dependent depth at q = 2fc Fj „ 
and q = fc F . „ + fc Pj „/ , respectively. The third contribution 
(Eq. A5) arises due to enforcing the constraint (25) on 
allowed deformations. 

The contribution to Eq. (A4) from evanescent modes 
with E n i > E F gives rise to the leading-order q- 
dependence of a non d, which is quadratic. This term 
therefore essentially captures the change of surface and 
curvature energy [cf. Eq. (11)]. Figure 5 compares the 
mode stiffness a, computed from Eqs. (A3)-(A5), to the 
Weyl approximation [first two terms on the r.h.s. of Eq. 



(11)]. Here R Q = h.7bk F 1 and L = lOOOfc/. Both 
the results for Au (i] s = 0.76, r\ c = —0.11) and for 
a pure constant-volume constraint (rj s = rj c = 0) are 
shown. In both cases, the overall g 2 -dependence of a 
for large q is evident, and consequently, the minimum at 
q = 2k F V :=y 0.56fcp is deeper than the other minima. 

When examining the softening of phonon modes, the 
overall increase of a{q) with q allows us to concentrate 
on the go = 2fc Fj v mode, where v labels the highest open 
channel. This mode will always be the dominant one. 
In general, the mode stiffness also shows less-pronounced 
minima at other values q > qo, due to the lower- lying 
channels, but these contributions are less singular than 
the go mode. It is therefore possible to approximate a 
by a smooth contribution tx q 2 plus the explicit terms of 
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0.8 
0.6 



<% 0.4 
H 



— r| s =0, 1 C =0 (const. 
r| s =0.76, Ti c =-0. 




FIG. 5: (Color online) The mode stiffness a computed from 
Eqs. (A3)-(A5) [solid curves]. The dashed curves show the 
Weyl approximation to a [first two terms on the r.h.s. of Eq. 
(11)]. The results for gold and for an idealized free-electron 
metal are shown; in each case, the radius is Ro = 5.75kp X and 
the length is L = lOOOfc^ 1 . 



Eq. (A3) for the highest open channel. 



APPENDIX B: PERTURBATION POTENTIAL 
MATRIX ELEMENTS 

Consider free electrons confined within a cylindrical 
nanowire of radius i?o and length L by a step poten- 
tial V{r) = VoO(r — Ro), where 9 is the step function. 
The transverse eigenf unctions in polar coordinates are 
given by if>^ n (r,ip) = (2tt)' 
dial function Xmnif) reads 



" 2 e mip Xmn{r), where the ra- 



J m ( V E mn r) , r < Rq, 

= N mn { J m (VE^;R )K m (^V„-E mn r) 
kJvV -E,~Ro) 



, r > Ro. 



(Bl) 



Here N mn is a normalization factor, J m is the Bessel 
function of order m, K m is the modified Bessel function 
of the second kind of order m, and for simplicity of nota- 
tion, we use the convention ft 2 /2m e = 1. The transverse 
eigenenergies E mn are determined by the continuity of 
drXmn/Xmn at r = Rq. An axisymmetric perturbation 
of the wire changes the confinement potential by 

8V(r, z) = V [6(r -R - 5R(z)) - 6{r - R )] , (B2) 

where the variation in radius is given by SR(z)/Rq = 
Eg t' q e lqz and the perturbation wavevectors q are re- 
stricted to integer multiples of 2ir/L. We expand the 
matrix elements of SV with respect to the unperturbed 



eigenfunctions * mn fe(r, (f, z) 
first order in b q and get 



L-VV^Jr,^) to 



(rank] SV\ mnk) 



V Ro (Ro)Xfhn(Ro) 

(B3) 



Taking the limit Vq — > oo we get 



El — e »(*-*+g)i 
q i(k 

q 



k + q)L 1 



lim (mnk\ SV\ mnk) ~ 2y E mn E„ 



% , 1 _ p i(k-k+q)L 



i(k — k+q)L 



In the limit of L — > oo we have 



lim (mnk\ SV\ mnk) ~ -I^J E mn E m n b q 5j, k _ q (B5) 



and recover a coupling of states with k = k + q, only. 
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